Topological and Transport Properties of Dirac Fermions in Antiferromagnetic 
Metallic Phase of Iron-Based Superconductors 
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We investigate Dirac fermions in the antifferomagnetic metallic state of iron-based superconduc- 
tors. Deriving an effective Hamiltonian for Dirac fermions, we reveal that there exist two Dirac 
cones carrying the same chirality, contrary to graphene, compensated by a Fermi surface with a 
quadratic energy dispersion as a consequence of a non-trivial topological property inherent in the 
band structure. We also find that the presence of the Dirac fermions gives the difference of sign- 
change temperatures between the Hall coefficient and the thermopower. This is consistent with 
available experimental data. 
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The iron-pnictide high-temperature superconductors 
have attracted much attention since the discovery 
In this multi-orbital system, the Fermi surfaces con- 
sist of electron and hole Fermi surfaces of compara- 
ble sizes. This multi-band structure leads to the rich 
physics, presumably including superconductivity with 
high-transition temperatures [H, [3 ■ In the parent antifer- 
romagnetic phase, the system remains metallic 0, Q con- 
trary to the simple antiferromagnetic state in single band 
systems. This parent antiferromagnetic state possesses a 
non-trivial topological property; consequently appears a 
Dirac fermion energy spectrum close to the Fermi energy 
0. In this Letter, we explore the Dirac fermion physics 
in this system. 

Dirac fermions in condensed matter systems were first 
highlighted by the observation of the anomalous integer 
quantum Hall effect in graphene 0, In graphene, the 
energy spectrum is linear at the corners of the first Bril- 
louin zone, so that the electrons at low energies are well 
described by the Dirac equation with the speed of light 
being replaced by the Fermi velocity 0]. Such a unique 
energy band structuregives rise to several characteristic 
transport properties [lO[. 

In the iron-pnictide superconductors, a similar linear 
spectrum was discussed by Ran et al. They pointed out 
that hybridizations between the Fe Sd orbitals and the 
pnictide ion 4p orbitals give rise to the band degeneracy 
characterized by a non-trivial topology. Contrary to a 
conventional spin density wave (SDW), there are gapless 
nodal points along the Fermi surface. The topology here 
is characterized by "vorticity" quantum number, which 
is associated with the phase winding defined through a 
two-component spinor wave function. The Fermi surfaces 
connected by the SDW wave vector have a vorticity zero 
and a vorticity two: This vorticity mismatch leads to 
a nodal SDW, and creates Dirac cones near the Fermi 



energy. 

In this Letter we shall derive the effective theory de- 
scribing the Dirac fermions based on the five-band model 
3 with the SDW mean field analysis 0, |Tl| . We reveal 
that there exist two Dirac cones carrying the same chi- 
rality contrary to the Dirac fermions in graphene. This 
unusual feature is due to the presence of another Fermi 
surface, which does not have a linear spectrum but has 
chirality. We demonstrate that the scattering rate differ- 
ence between Dirac fermions and conventional electrons 
leads to anomalous temperature dependence of the trans- 
port coefficients which are consistent with the experiment 

For the band structure calculation, we take the five 
band model 
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where dj^^, creates an electron with spin s on the fj,-th 
orbital at site j. The parameters ti^jv and are given 
in Ref . 3 . Following Ref . [1] we take the following inter- 
action form: 

-Hi = Uj^^it^t^Jt^i + i^-'^J) ^msrijus' 

j,fj,<u,s,s' 

'^It.sd'lus'djf^s'djus 

j,fj,<iy.,s,s' 
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where n,„s = dl,,d\,^^. The mean field Hamiltonian for 

J J ^S J ^s 

the SDW state reads, 



'Electronic address; |morinari@yukawa. kyoto-u.ac.jp] 



l<i£RBZ,fi,u,s 



d^ 



[n 



ks 



fiiy 



(3) 



2 




• Fe 'As OCa _0.5eV 



FIG. 1: (color online) (Left) The crystal structure of 
CaFe2As2. (Middle) The electron band structure near the 
Fermi energy in the SDW state. The unit cell is taken by the 
square lattice formed by Fe atoms. (Right) The magnification 
of the band structure near the Dirac point indicated by the 
circle in the left. 

where Q = (tt, 0) and the k summation is taken over the 
redueed Brillouin zone. The matrix "Hies is given by 

["^ks]^,, = ,Jfiy tS°^ ) , (4) 

where we take the plus (minus) sign for s =t (i). The 
mean field order parameters, m^j^-i^, are determined by 
solving the mean field equations in a self-consistent way. 
Figure [T] shows the electron band structure calculated in 
the SDW state with f/ = 1.2 eV and J = 0.25 eV with the 
zero temperature magnetic moment 1.0/iB, which is close 
to the value observed by neutron scattering in BaFe2As2 

In the SDW state, Dirac points appear near the Fermi 
energy. Those Dirac points exist on the kx axis and, 
in some parameter region, on the ky axis as well but 
not at the high symmetric points in the Brillouin zone. 
Since the Dirac points appearing on the ky axis are not 
protected by topology @, we focus on the Dirac points 
on the kx axis. These Dirac points are stable as long as 
the magnetic moment is less than i^B- 

In order to investigate the physical properties of 
the Dirac fermions we derive the effective Hamiltonian 
around the Dirac point following Ref. p^ . in which 
the Dirac fermions in the organic conductor a-(BEDT- 
TTF)2l3 were analyzed. The effective Hamiltonian is 
used to determine chirality and the Landau level struc- 
ture of Dirac fermions. 

First we expand the Hamiltonian "Hk around the Dirac 
point kc with denoting k = kc -I- <5k, 

= Hk. + ^-^5kx + ^6ky + ... (5) 

(The spin index is suppressed hereafter.) Next, noting 
that a Dirac point is formed by two eigenstates we con- 
struct the effective 2x2 Hamiltonian using those two 
eigenstates denoted by |k, ±) . The general form of the 
effective Hamiltonian is [IJ] , 

n'sU= E '^k-v.ap, (6) 



where a = {a^ ,ay,az) are the Pauh matrices and ctq is 
the unit matrix. (We set H = 1.) The parameters Vp are 
determined by the following equation with a = x,y: 

p—Q,x,y,z 

Figure [2][a) shows the energy band dispersion of the 
state (shown in Fig. [1]) along the k^ axis. The Dirac 
points are at kc = (±0.829, 0). The positions of the Dirac 
points are close to the values reported in angle resolved 
photoemission spectroscopy measurements [l5|. The ef- 
fective Hamiltonian is given by 

= {Ec ± v^kx) (To ± {v^^kxcr, + vlkyUx) , (8) 

where - Ep ^ ~20 meV, v^/a = 0.286 eV, vl/a = 
0.229 eV, /a = 0.672 eV with a the lattice constant. 
The non-zero value of Ug implies that the Dirac cone is 
tilted in the k^ direction. As a result, we expect a strong 
anisotropy in the Fermi velocity. 

Now we discuss chirality associated with Dirac 
fermions. In Fig. [DJb) and (d), the vector nk = 
(k, -|-|<T|k, is shown around the two Dirac points. The 
vector field nk shows a vortex configuration with vor- 
ticity one. A remarkable fact is that the two Dirac 
cones have the same vorticity: The vector nk is ro- 
tated clockwise when one goes around each Dirac point. 
In conventional Dirac fermion systems, like graphene, 
two Dirac cones have the opposite chirality by symme- 
try. The same is true for the Dirac cones in a-(BEDT- 
TTF)2l3. This unusual property is understood by consid- 
ering the chirality associated with the hole band around 
r point as shown in Fig. ^c). Although the energy 
dispersion around the F point is quadratic, and an en- 
ergy gap lies between two relevant bands, we can apply 
the above analysis in constructing the effective Hamilto- 
nian: In the quadratic dispersion case, the parameters 
Vp are linear functions with respect to k^ or ky with 
higher order corrections. Along the circle represented by 

k = ^Jk^~+lt^ = 0.1, we find V.^ = flk^o + b^a^ + c^ctx, 
with flk = -0.508fc2 - 0.475fc2, 6k = -0.191fc2 + 2.56fc2, 
and Ck = 2.80kxky. The winding is described by the fol- 
lowing vector: {nx,nz) = {ck,bk) / \/l>i, + c^. As shown 
in Fig.[2][c), the vector nk rotates counter-clockwise twice 
when one goes around the F point. This chirality exactly 
compensates chirality of the two Dirac cones. We note 
that the scattering between the two Dirac cones, which 
is still under debate in graphene [l6l | , is qualitatively dif- 
ferent from that in graphene because these Dirac cones 
carry the same chirality. It is worth of pointing out that 
this chirality property is compatible with the interband 
Cooper pairing between the hole band around the F point 
and the Dirac cones. 

Now we turn to the physical properties associated with 
the Dirac fermions. One of the most significant effects as- 
sociated with the Dirac fermion chirality is the suppres- 
sion of backward scattering [l3|. As a consequence the 
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FIG. 2: (color online) (a) The energy bands along the kx 
axis. The dotted line represents the Fermi energy. Two Dirac 
points are located at (±0.829, 0). The direction of the vector 
n(k) = (n^(k),n^(k)) around (a) (-0.829,0), (b) (0,0), and 
(c) (0.829,0). 

scattering ratio is different between Dirac fermions and 
conventional electrons. Reflecting this difference some 
transport coefficients clearly exhibit contributions from 
Dirac fermions If 811 even if the Dirac fermions are the 
minor carrier [l9|. 

We consider a phenomenological two-band model con- 
sisting of a hole band, with a conventional energy spec- 
trum, and an electron band, with the Dirac fermion 
energy spectrum. The charge current is given by j = 
(tE -I- a(— VT), where a and a are tensors computed 
by the Jones- Zener solution of the Boltzmann equation 
[20|. The contribution from the Dirac fermions is de- 
noted by cr'-'^-' and a*-"^-* and that from the holes is de- 
noted by a^^^ and a^'^\ The sign change of the Hall 
coefficient Rh = '^xy/iBzc'^x) with the applied mag- 
netic field in the z direction occurs at aiij /ai'y' = 1. On 
the other hand, the sign change in the thermopower oc- 
curs when (x'xl joL^xx = 1- For the energy dispersions, we 
take e]^ = vk — Se and e^'*-' = —k'^/(2m) + eh with energy 
levels Ee and eh . At finite temperatures, electrons may be 
scattered by phonons and spin waves. Here we focus on 
temperatures lower than the characteristic temperatures 
associated with those excitations and assume constant 
scattering times, Te.h- In order to reduce the number of 
parameters, we set mv'^/2 = Be = Sq. The difference 
between Tg and is parameterized by = t^It^. In 

Fig. [3]^a) and (b) we show ^/rjr^ = \fii^^^J{T^a^) 

and Se/sh = T^a^xx / {teCx^J) for differnt values of rir = 
Tie/ Ufi with n^(^h) the electron (hole) density. Whether 
there are sign changes or not depends on and n^. (Note 
that Ur is controlled by changing e/j.) The sign change 
occurs when = r^/rh, Se/sh- In general, the sign 
change occurs at different temperatures for the Hall co- 
efficient and the thermopower because of their different 
dependences on r^. 



We show the Seebeck coefficient, the Hall coeffi- 
cient, and the Nernst coefficient, v = {uxx^^xy ~ 
^xy^xx )/ [Bz{(^Ix + '^ly)] 7 in Fig- M.c) for = 0.45 and 
Ur = 0.05. The result is consistent with the experiment 
reported in Ref . [l2j , where the similar sign changes were 
observed in CaFe2As2. Note that we expect large contri- 
butions from the Dirac fermions for the Nernst coefficient 
as well because the quantity is quadratically dependent 
on the scattering time. We point out that a naive appli- 
cation of the Jones-Zener solution to the SDW mean field 
state with the same scattering time for all bands does not 
lead to these temperature dependence even if we take into 
account the temperature dependence of the order param- 
eter. In CaFe2As2, the energy dispersion in the kz direc- 
tion may not be negligible. Therefore, we assumed here 
a moderate value for t^. For systems with strong two- 
dimensionality, like the so-called 1111 system, we expect 
a small t^. In such a case Dirac fermions markedly con- 
tribute to the transport coefficients. 

If we assume a conventional energy dispersion f or the 
electron band, the temperature dependence of -^/r^Jrh 
is qualitatively different from the Dirac fermion case as 
shown in the inset in Fig. [2Ia). In this case, the Hall 
coefficient behaves qualitatively different from Fig. |31Jc). 
The Hall coefficient is positive at low temperatures, unlike 
the Dirac fermion case. 

A systematic change of these sign changes are expected 
by controlling the electron density. In Fig.|4]we show the 
doping dependence of the Dirac point and its position in 
the Brillouin zone in the five-band SDW state. Both de- 
crease monotonically by increasing electron density. In- 
terestingly the Dirac point crosses the Fermi energy at 
the hole doping side. 

Now we propose an experiment to directly observe the 
evidence of the Dirac fermions. In conventional electrons, 
the Landau level energies are equally separated. By con- 
trast, the Dirac fermion Landau level energies are not. 
For the tilted Dirac cone (|8]), the Landau level energy 
spectrum is given by [2l| 



E,, = sgn{n)^^^^/2)^\, (9) 

with £b ~ XjsfeB the magnetic length and A = 
-\/l — (wq/uJ)2. This Dirac fermion Landau level struc- 
ture can be observed by scanning tunneling spectroscopy 
(STS) measurements. At graphite surfaces under mag- 
netic fields a series of peaks in the tunnel spectra as- 
sociated with Landau quantization of the quasi-two- 
dimensional electrons were observed [1^ [2^. The ob- 
served Landau levels are in good agreement with those 
expected for a surface layer, graphene. From our effec- 
tive Hamiltonian, we evaluated that E\ ~ 10 meV at 
S = 10 T while Ei ~ 130 meV in graphene. Using STS, 
we expect that the Dirac fermion Landau level structure 
should be observed in the parent antiferromagnetic state. 

To conclude, we have investigated the Dirac fermions 
in the parent antiferromagnetic state of the iron-based 
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superconductors. There exist two tilted Dirac cones 
with the same chirahty contrary to the Dirac cones in 
graphene. The sign changes observed experimentally in 
the thermopower and the Hall coefRcient are consistently 
understood from the Dirac fcrmion picture. We propose 
that the Dirac fcrmion Landau level structure be ob- 
served by STS measurements. Since the two Dirac cones 
carry the same chirality, the intervalley scattering is qual- 
itatively different from ordinary Dirac fermion systems. 
Compared to graphene, the valley splitting occurs in a 
different manner that would be clarified in the magne- 
toresistance. 
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FIG. 4: (color online) The electron density dependence of 
the Dirac point and the Dirac point position in the Brillouin 
zone denoted by (fca;, 0). The inset shows the electron density 
dependence of the magnetic moment. 
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FIG. 3: (color online) The temperature dependence of the ra- 
tios ^re/rh (a) and Se/sh (b) for different values of the den- 
sity ratio nr. {ur = f.O, 0.70, 0.45, 0.35, 0.25, O.f 5, O.fO, 0.05 
from the top to the bottom.) The inset in (a) shows the 
temperature dependence of yrv/rVt for conventional energy 
dispersions, {rir — f.O, 0.45, 0.25, O.f 0, 0.05 from the top to 
the bottom.) (c) The temperature dependence of the Hall 
coefficient, the Seebeck coefficient, and the Nernst coeffi- 
cient, which are normalized by (27ra^/|e|)(to/eofl)'^, ks/le], 
and {kBO? / {2-Kh)){hv /£oa){vTD / a), respectively, with fcs the 
Boltzmann constant. 
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